addpath('npy-matlab-master\npy-matlab')
A2 = readNPY('CPGE_map_680nm_-24V_30K_y.npy')
A2=A2*1e12
for i =1:16
    for j =1:16
        if abs(A2(i,j))>2
            A2(i,j)= (A2(i,j-1)+A2(i,j+1))/2
        end
    end
end
A2=fliplr(fliplr(A2).')
surf(A2,'edgecolor','none')
set(gca,'xlim',[6 16],'ylim',[2 12])
view(2)
%colormap(winter)
caxis([-2.2,2.2])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
colorbar('Ticks',[-2:2:2],'TickLabels',{})
set(gcf,'Position',[0, 0, 380, 300])